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The scalable application of quantum information science will stand on 
reproducible and controllable high-coherence quantum bits (qubits). Here, we 
revisit the design and fabrication of the superconducting flux qubit, achieving a 
planar device with broad frequency tunability, strong anharmonicity, high 

reproducibility, and relaxation times in excess of 40 ps at its flux-insensitive point. 
Qubit relaxation times T x across 22 qubits are consistently matched with a single 

model involving resonator loss, ohmic charge noise, and 1/f flux noise, a noise 
source previously considered primarily in the context of dephasing. We 
furthermore demonstrate that qubit dephasing at the flux-insensitive point is 
dominated by residual thermal photons in the readout resonator. The resulting 
photon shot noise is mitigated using a dynamical decoupling protocol, resulting in 
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T 2 « 85 jus, approximately the 27J limit. In addition to realizing an improved flux 
qubit, our results uniquely identify photon shot noise as limiting T 2 in 

contemporary qubits based on transverse qubit-resonator interaction. 


Over the past 15 years, superconducting qubits have achieved a remarkable five-order- 
of-magnitude increase in their fundamental coherence metrics, including the energy 
decay time T\ , the Ramsey free-induction decay time 73* , and the refocused Hahn- 
echo decay time 73 e- This spectacular trajectory is traceable to two general strategies 
that improve perfonnance: (1) reducing the level of noise in the qubit environment 
through materials and fabrication improvements, and (2) reducing the qubit sensitivity 
to that noise through design advancements 1 . 

The charge qubit evolution is a quintessential example 2 . Early demonstrations (Cooper- 
pair box) exhibited nanosecond-scale coherence times 3 . Since then, operation at noise- 
insensitive bias points (quantronium) 4 , the introduction of capacitive shunting 
(transmon) 5 , the use of 2D 6 and 3D 7 resonators to modify the qubit electromagnetic 
environment, the development of high-Q capacitor materials and fabrication 
techniques 8,9 , and the introduction of alternative capacitor geometries (Xmon) 10 have 
incrementally and collectively raised coherence times to the 10-100 ps range 10,11 and 
beyond 12,13 . In addition, the capacitive shunt has generally improved device-to-device 
reproducibility. The trade-off, however, is a significant reduction in the charge qubit 
intrinsic anharmonicity (i.e., the difference in transition frequencies / 01 and f n between 
qubit states 0, 1 and 1, 2) to 200-300 MHz for contemporary transmons, complicating 
high-fidelity control and exacerbating frequency crowding in multi-qubit systems 14 . 

In contrast, the perfonnance of the persistent-current flux qubit 15,16 has progressed more 
slowly over the past decade. Device asymmetry was identified early on to limit flux 
qubit coherence 17 and, since 2005, symmetric designs have generally achieved 0.5-5 
ps 18,19 with a singular report of T 2E - 23 jus ~ 25j 20 . Despite respectable performance 

for individual flux qubits, however, device-to-device reproducibility has remained poor. 
An early attempt at capacitive shunting 21 improved reproducibility, but coherence 
remained limited to 1-6 //s 22,23 . Recently, flux qubits embedded in 3D 24 and coplanar 25 
resonators exhibited more reproducible and generally improved relaxation and 
coherence times: 7j = 6 - 20 jus , 7V =2-8 //s . Nonetheless, further improvements in 

these times and in reproducibility are necessary if the flux qubit is to be a competitive 
option for quantum information applications. 

In this context we revisit the design and fabrication of the flux qubit. Our 
implementation, a capacitively shunted (C-shunt) flux qubit 21 coupled capacitively to a 
planar transmission-line resonator, exhibits significantly enhanced coherence and 
reproducibility while retaining an anhannonicity varying from 500 - 910 MHz in the 
four devices with the highest relaxation times. We present a systematic study of 22 
qubits of widely varying design parameters - shunt capacitances C sh = 9 - 51 fF and 
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circulating currents I =44- 275 nA - with lifetimes at the flux-insensitive bias point 
ranging from 7j < 1 /us (small C sh , large I ) to T x = 55 jus (large C sh , small I ). Over 
this entire range, the measured 7j values are consistent with a single model comprising 
ohmic charge noise, 1/ f flux noise, and Purcell-enhanced emission into the readout 
resonator. We furthermore investigated and identified quasiparticles as a likely source 
of observed T x temporal variation. For the highest coherence devices, the Hahn-echo 

decay time T 2E = 40//s < 2T X does not reach the 2T X limit, as is also often observed with 
transmons coupled transversally to resonators 710,26 . We demonstrate that this is due to 
dephasing caused by the shot noise of residual photons in the resonator (mean photon 
number n 0 = 0.006), observing a lorentzian noise spectrum with a cut-off frequency 
consistent with the resonator decay rate. We then use Carr-Purcell-Meiboom-Gill 
(CPMG) dynamical decoupling to recover 7) C|)MG « 2 T x in a manner consistent with the 
measured noise spectrum. 

Results 

C-shunt flux qubit. Our circuits each contain two C-shunt flux qubits—with different 
frequencies—placed at opposite ends of a half-wavelength superconducting coplanar 
waveguide (CPW) resonator (Fig. la). The resonator, ground plane, and capacitors (Fig. 
lb) were patterned from MBE-grown aluminum deposited on an annealed sapphire 
substrate 8 (Supplementary Note 1). We used both square capacitors (Fig. lb) and 
interdigital capacitors (IDCs, not shown) coupled capacitively to the center trace of the 
CPW resonator to enable qubit control and readout. In a second fabrication step, the 
qubit loop and its three Josephson junctions (Fig. lc) were deposited using double¬ 
angle, electron-beam, shadow evaporation of aluminum. One junction is smaller in area 
(critical current) by a factor a , and each of its leads contacts one electrode of the shunt 
capacitor. An equivalent circuit is illustrated in Fig. Id (Supplementary Note 2). 

Varying the qubit design enables us to explore a range of qubit susceptibilities to flux 
and charge noise with impact on both T x and T 2 21 . Compared with the conventional 

persistent-current flux qubit 15,16 , our best C-shunt flux qubits have two key design 
enhancements. First, a smaller circulating current - achieved by reducing the area and 
critical current density of the Josephson junctions (Fig. lc) - reduces the qubit 
sensitivity to flux noise, a dominant source of decoherence in flux qubits. Second, a 
larger effective junction capacitance - achieved by capacitively shunting the small 
junction (Fig. lb) - reduces the qubit sensitivity to charge noise, and improves device 
reproducibility by reducing the impact of both junction fabrication variation and 
unwanted stray capacitance. Furthermore, the use of high-quality fabrication techniques 
and physically large shunt capacitors reduces the density and electric participation of 
defects at the various metal and substrate interfaces 1 . 

The system is operated in the dispersive regime of circuit quantum electrodynamics 
(cQED) and is described by the approximate Hamiltonian 27 

n ilv «h % (<S> i )& I l2+h t o,(n + m)+2h X (t> t )(n + \/2)&j2, ( 1 ) 
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where the three terms are respectively the qubit (represented as a two-level system), 
resonator, and qubit-resonator interaction Hamiltonians, <r z is the Pauli operator 

defined by the qubit energy eigenbasis, co r is the resonator angular frequency and h is 

the resonator photon-number operator. The qubit angular frequency &> q (® b ) is set by 

the magnetic flux bias O b , measured relative to an applied flux (m + 1/2)0 0 where m 

is an integer and O 0 is the superconducting flux quantum, and attains its minimum 

value 0) q (O) = A at the flux-insensitive point O b = 0. The quantity/(O b ) is the qubit- 

state-dependent dispersive shift of the resonator frequency, which is used for qubit 
readout. In the Supplementary Notes 3-5, we discuss further the two-level-system 
approximation for the C-shunt flux qubit, an approximate analytic treatment which goes 
beyond Eq.(l), and the numerical simulation of the full qubit-resonator Hamiltonian 
used to make quantitative comparisons to our data. 

7j relaxation and noise modeling. We begin by presenting the T ] characterization 
protocol for the device in Fig. 1. We first identify the resonator transmission spectrum 
(Fig. 2a, top panel) by scanning the readout-pulse frequency co m about the bare 

resonator frequency cojIn ~ 8.27 GHz. Using standard cQED readout, qubit-state 
discrimination is achieved by monitoring the qubit-state-dependent transmission 
through the resonator 27 . Next, we add a qubit driving pulse of sufficient duration to 
saturate the ground-to-excited-state transition and sweep the pulse frequency (O d (Fig. 
2a, bottom panel). The resulting spectra for qubits A and B (Fig. la) exhibit minima 
A a /2 n ~ 4.4 GHz and A B /2 n « 4.7 GHz at the qubit flux-insensitive points and 
increase with magnetic flux (bias current) away from these points. Finally, using a 
single n -pulse to invert the qubit population, we measure the 7j relaxation of qubit A ( 

T ] = 44//s ) and qubit B (7j = 55//s ) at their flux-insensitive points (Fig. 2b). High- 

power spectroscopy (see Supplementary Note 6) reveals transitions amongst the first 
four qubit energy levels that are well matched by simulation, and identifies 
anharmonicities of 500 MHz in the two measured devices. 


Using this protocol, we investigated 22 C-shunt flux qubits from five wafers 
(fabrication runs), spanning a range of capacitance values (C sh = 9 - 51 fF ) and qubit 
persistent currents (/ = 44-275 nA) and featuring two capacitor geometries 
(interdigital and square). The junction critical currents were adjusted to maintain 
A/2;r « 0.5 -5 GHz (see Supplementary Note 7). 


The data were analyzed using simulations of the full system Hamiltonian and a Fermi’s 
golden rule expression for the exited state decay rate 21 , 
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where | g'j (| e) ) indicates the qubit ground (excited) states, and the sum is over four 

decay mechanisms: flux noise in the qubit loop, charge noise on the superconducting 
islands, Purcell-enhanced emission to the resonator mode, and inelastic quasiparticle 

tunneling through each of the three junctions. The operator D x is a transition dipole 
moment, and 5/ ( co l] j is the symmetrized noise power spectral density for the 

fluctuations which couple to it. For example, D (|1 is a loop current operator for flux 
noise (&>) , and D 0 is an island voltage operator for charge noise S Q (&>) 
(Supplementary Note 8 and 9). 

We considered both S- (co) qc 1 / of (inverse-frequency noise) and S x ( co ) oc co (ohmic 

noise) - the two archetypal functional forms of noise in superconducting qubits 20,28 " 33 - 
for our magnetic flux and charge noise models, and used the frequency dependence of 
7j for specifically designed devices to distinguish between them. While the following 

results are presented using symmetrized power spectral densities, we are careful to 
account for the distinction between classical and quantum noise processes in making 
this presentation (Supplementary Note 9). 

For example, in Fig. 3a, Qubit C (C sh = 9 fF) has a large persistent current 
(/ =275 nA) and a small qubit frequency ( A c /2;r = 0.82 GHz), making it highly 
sensitive to flux noise. Consequently, the measured 7j is predominantly limited by flux 
noise over a wide frequency range. This 7)-trend constrains the flux noise model to the 
form Sq (ry) = A 2 (2^x1 Wz/co) 7 over the range 0.82 - 3 GHz (black dashed line, Fig. 
3a). For comparison, the functional form for ohmic flux noise (grey dashed line), scaled 
to match 7j at A c /2;r = 0.82 GHz (green dot), is clearly inconsistent with all other data 

over this frequency range. The noise parameters A 2 = (l .4 /kl> 0 )~ /Hz and y = 0.9 used 

to match the data in Fig. 3 a are derived from independent measurements - Ramsey 
interferometry 31 and T\ p noise spectroscopy 32 (Supplementary Note 10) - made at much 
lower frequencies in the context of classical noise related to qubit dephasing (Fig 3b). 
These values are commensurate with earlier work on qubits 20,31 ' 33 and dc 
Superconducting QUantum interference Devices (SQUiDs) 34 . The consistency between 
the magnitude and slope of the flux noise power spectra, spanning more than twelve 
decades in frequency - millihertz to gigahertz - is remarkable, made even more so by 
the fact that the data in Fig. 3b were measured with a different device (qubit B, Fig.3c). 

In contrast, Qubit B (C sh = 51 fF) has a much smaller persistent current (I =49 nA ) 
and larger qubit frequency ( A |s / 27T = 4.7 GHz ). Its value of JJ is most strongly 
influenced by charge noise (magenta dashed line, Fig. 3c) in the 5.0 - 6.5 GHz range, 
consistent with an ohmic charge noise model of the fonn S 0 (&>) = 4 2 coj (in x 1 GHz) 

with the parameter = (5.2x10 9 e) 2 /Hz adjusted to match the data. In addition to flux 
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and charge noise, the predicted value of 7j due to Purcell loss (light blue dashed line) is 
also included in Fig. 3a and 3c and involves no free parameters (see Supplementary 
Note 8). The resulting net value of T ] due to all three mechanisms (inverse-frequency 
flux noise, ohmic charge noise, and Purcell loss) is indicated with a red solid line and is 
in relatively good agreement with the ceiling of measured T { values. As we describe 
below, quasiparticles are responsible for reducing the 7j below this ceiling. 

Using these models, Fig. 3d shows a comparison of the measured and predicted T { 
values for all 22 qubits. The flux noise model (from Figs. 3a and 3b) is applied to all 
qubits, and the Purcell loss is included with no free parameters. For the charge noise 
model, to achieve agreement across all devices, it was necessary to use 
Aq sq =(5.2x1 0 9 c) 2 /Hz for square capacitors (from Fig. 3b) and 

Alq idc = (11.0 x 10 9 c) 2 /Hz for interdigital capacitors (IDCs), presumably reflecting the 

larger electric participation of the surface and interface defects for the IDC geometry 1 . 
The agreement is noteworthy, given that these qubits span a wide range of designs 
across five fabrication runs (see Supplementary Note 7). 

We note that inverse-frequency charge noise was incompatible with these data over the 
entire frequency range investigated (not shown), implying that the cross-over between 
inverse-frequency and ohmic charge noise occurred at a frequency below 0.82 GHz. 
However, while ohmic flux noise S# (co) cc co was inconsistent with T x over the 

frequency range 0.82 - 3 GHz, its functional form is plausibly consistent with data 
above 3 GHz when appropriately scaled (upper dashed grey line, Fig. 3a) and, therefore, 
cannot be conclusively distinguished from ohmic charge noise. Although the best 
agreement across all 22 qubits (Fig. 3d) did not require ohmic flux noise, we could not 
rule out its presence in the 3-7 GHz range. In Supplementary Note 11, we compare 
models that use ohmic charge noise (as in Fig. 3) and ohmic flux noise. Differentiating 
between such charge and flux noise at higher frequencies will be the subject of future 
work. Indeed, for both ohmic flux noise S m (co) cc co and inverse-frequency charge 

noise S Q (co) oc \jco , it is certainly possible (even expected) that the former (latter) 
dominates the flux (charge) noise at sufficiently higher (lower) frequencies. 

The measured data for qubit B (Fig. 3c) exhibit fluctuations in the range 7j = 20 - 60//s 
for qubit frequencies co^jln = 4.7 - 6.5 GHz . To investigate their temporal nature, we 
measured 7j repeatedly at the qubit flux-insensitive point co^jln = A B /2;r = 4.7 GHz 

over a 10-hour period and collected the data into sets of 50 individual decay traces. 
Figures 4a and 4b show the results of two such experiments, with set 2 being taken 
approximately ~ 17 h after set 1. The average of all traces from set 1 exhibits a purely 
exponential decay, whereas the corresponding average for set 2 exhibits a faster short- 
time decay and clear non-exponential behavior (Fig. 4a). Histograms of the T x values 
for individual traces exhibit a tight, Gaussian-shaped distribution centered at 55//s for 
set 1 and a broader, quasi-uniform distribution centered near 45//s for set 2. Over the 
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course of several weeks, we observed transitions between these two characteristic 
modes of behavior every few days for this device 35 . 

We attribute the temporal fluctuations and non-exponential decay function to excess 
quasiparticles - above the thennal equilibrium distribution - near the qubit junctions 36 " 

39 . Following Ref. 40, we define 7] as the average relaxation time associated with a 

single quasiparticle and take the quasiparticle number n qp to be Poisson-distributed with 
mean value n qp . This results in a qubit polarization decay function, 


(P e (0) = e^ (exp( -^ ) - l) e-^ R , (3) 

where 7j R captures the residual exponential decay time in the absence of quasiparticles 
(n = 0 ). The non-exponential decay function observed for set 2 is well described by 

Eq. (3) [black line in Fig. 4a] with fitting parameters n qp = 0.26 , T lqp = 23//S and 
T lR - 60 /js . 

We use a quantum treatment of quasiparticle tunneling to model the impact of single 
quasiparticles on the 7] of qubit B (Supplementary Note 8). Using a quasiparticle 

density x qp =4x10 (per superconducting electron), the calculated T lqp recovers the 

fitted value T lqp = 23 //s at the flux-insensitive point. Both T lqp and x qp are comparable 

to the quasiparticle-induced relaxation rates and quasiparticle density reported for 
similar devices 24 41 . The shaded region in Figs. 3a and 3c indicates the range of 
predicted 7] in the presence of n qp =0-1.0 quasiparticle. Most 7] data he within this 
region, supporting the hypothesis that their scatter (particularly for qubit B in Fig. 3c) 
and the observed temporal 7] variation (Fig. 4b) arise from the common mechanism of 

quasiparticle tunneling. Additionally, the residual relaxation time 7[ R for set 2 is similar 
to the exponential time constant obtained for set 1, indicating an underlying consistency 
in the noise models between the two data sets in the absence of quasiparticles. Unlike 
qubit B, qubit C consistently exhibited an exponential decay function (Fig. 4c) with 
little temporal variation (Fig. 3a, 4d), indicating that quasiparticles did not strongly 
influence this device. 

The results of Fig. 3 and Fig. 4 demonstrate clearly that 1/ f -type flux noise is the 
dominant source of qubit relaxation for frequencies below 3 GHz. To further strengthen 
this claim, it is instructive to compare relaxation times for qubits with similar 
frequencies and shunting capacitances, but where the persistent current (and thereby the 
sensitivity to flux noise) differs. We find that by reducing I from 170 nA to 60 nA, we 

improve the measured 7] from 2.3 to 12 us (see qubits 11 and 13 in Supplementary 
Table 1 in Supplementary Note 7). 
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Pure dephasing and thermal photon noise. We now address the transverse relaxation 
timer, and our ability to refocus coherent dephasing errors. Efficient refocusing 

implies that T 2 is limited entirely by 7], since 1/ T 2 = 1/27] +1 jT^, where T tp is the 
dephasing time. Generally, T 2 is maximal at the flux-insensitive point for conventional 
flux qubits 18 " 20 , and the device reported in Ref. 20 was efficiently refocused with a 
single echo pulse ( T 2E - 23 //s « 27]). In the current work, however, a single refocusing 
pulse is no longer completely efficient ( T 2E < 27]). This suggests that an additional, 

higher-frequency noise channel has been introduced. Unlike the device in Ref. [20], 
which was coupled to a dc SQUID for readout, our C-shunt flux qubits are transversally 
coupled to a resonator (Fig. 1). Such inefficient refocusing is also reported for 
transmons similarly coupled to resonators 710,26 . 

As we show below, the main source of dephasing in C-shunt flux qubits biased at their 
flux-insensitive point is photon-number fluctuations (shot noise) in the resonator, which 
vary the qubit frequency via the ac Stark effect (as in the transmon case 26,27 ). Given a 
small thermal-photon population n = 1 / [e tm ' lksT -1) <SC 1 in the resonator (see 
Supplementary Note 12), the photon-induced frequency shift A* ark and dephasing rate 
fj of the qubit are 42 


A stark='/ 2 Z« , 

(4) 

r t=v—n- 

K 

(5) 


The factor rj = tc 2 / (k + 4% 2 ) effectively scales the photon population seen by the 
qubit due to the interplay between the qubit-induced dispersive shift of the resonator 
frequency % and the resonator decay rate k . Both the strong dispersive (2 % » k) and 
weak dispersive {2% «: k) regimes have been previously addressed 26,43,44 . Here, we use 
qubit B to focus primarily on the intennediate dispersive regime {2x/2k = 0.9 MHz, 
at/ 2zr = 1.5 MHz, see Fig. 5a) relevant for high-fidelity qubit readout 45 . 

We begin by intentionally injecting additional thennal photons n M {P M ) into the 
resonator from an external noise generator with power P add (Fig. 5b and Supplementary 
Note 2). In the small- n add limit, the measured qubit spectrum exhibits a linear 

relationship between the effective qubit frequency (o' = &> q + A Lamb + A* ark (n add ) and the 
generator power P add (Figs. 5c and 5d). For completeness, we have included the Famb 

shift A Lamb , a fixed frequency offset due to the resonator zero-point energy. Combining 
the extracted slope with Eq. (4), we calibrate the dependence of the added-photon 
population n add (in the resonator) on the generator power 7] idd . 



Next, we measure the Hahn-echo dephasing rate for several photon populations using 
the calibrated n. M (P add ). All echo traces (Fig. 5e) feature exponential decay rates 

r 2E = \/T 2e , indicating little (if any) impact from 1// noise (charge, flux, ...) and 
consistent with photon shot noise featuring a short correlation timel/x" <SC T 2E . The 
extracted pure dephasing rate T E = T 2E -T)/2 scales linearly with photon population 
n add (P ad d) (Fig. 5f). The extracted slope agrees with Eq. (5) to within5%. The non-zero 
dephasing rate at n aAd - 0 corresponds to a residual photon population n 0 = 0.006, 
equivalent to an effective temperature T e{( = 80 mK. By comparison, the qubit effective 
temperature determined from its first excited-state population is 35 mK 13 . 

To confirm that the noise arises from residual thennal photons, we directly measure the 
noise power spectral density (PSD) using the T lp (spin-locking) method 32 . This method 

(inset Fig. 6a) collinearly drives the qubit along the Y-axis with a long Y pulse, which 
“locks” the qubit state in the rotating frame. Measuring the qubit relaxation rate in the 
rotating frame, r ij0 (Q Rabi ) = £ z (Q Rabi )/ 2 + T l /2 , effectively samples the noise PSD 

S z {co) seen by the qubit at the locking (Rabi) frequency 0 Rabi (see Supplementary Note 
10). By varying the locking drive amplitude, which is proportional to Q Rabi , we sample 
the noise spectrum over the range®/ 2k = 0.1 - 100 MHz (Fig. 6a). Below 10 MHz, the 
resolved noise spectra for all n add (including« add = 0) have similar shapes: flat (white) 
at low frequencies with a 3-dB high-frequency cut-off at the resonator decay rate co — k 
. This fonn is consistent with the expected lorenzian PSD for thennal-photons in a 
resonator as seen by the qubit (see Supplementary Note 10), 


S z (a) = (2 x f 


IrjnK 
co 2 +K 


2 ’ 


( 6 ) 


which includes the dispersive coupling % and the filtering factor // [see Eqs. (4,5)]. 
Equation (6) agrees with the measured PSDs for all photon populations n - n aAA + n 0 , 
with the residual photon number n 0 extracted from Eq. (6). This agreement eliminates 

the driving or readout field as the source of the residual photons, because such coherent- 
state photons follow Poisson statistics with a resulting cut-off frequency k / 2 (one-half 
the observed value) 46,47 . 

Finally, we apply dynamical decoupling techniques to validate the functional form of 
the measured noise PSD and to recover T 2 « 27]. We use the Carr-Purcell-Meiboom- 

Gill (CPMG, inset Fig. 6b) pulse sequence, comprising a number N n of equally spaced 

n -pulses. The application of n -pulses in the time domain can be viewed as a bandpass 
filter in the frequency domain which shapes the noise spectra seen by the qubit 21,48 " 50 . 
Since the filter passband is centered at a frequency inversely related to the temporal 
spacing At between adjacent pulses, increasing N n for a fixed sequence length will shift 
this passband to higher frequencies (see Supplementary Note 13). 
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Figure 6b shows the measured CPMG decay time T 2( mG vs. it -pulse number N K with 

no added noise (« add = 0 ). From N K = 1 (Hahn echo) to N K =100, the decay time 

^ 2 cpmg remains near 40 /us , consistent with the white-noise (flat) portion of the noise 

PSD in Fig. 6a. Above N n - 100, the passband frequency traverses the cut-off region of 

the PSD and, as the integrated noise level decreases, 7) CPMG rises. For N n > 1000, the 

refocussing becomes efficient with f 2CPMG « 85//s « 27j. The close correspondence 

between the noise spectral density in Fig. 6a and the mitigation of that noise by CPMG 
in Fig. 6b strongly supports our methods and interpretations. 

Discussion 

The C-shunt flux qubit is a planar device with broad frequency tunability, relatively 
strong anharmonicity and high reproducibility, making it well suited to both gate-based 
quantum computing and quantum annealing. The anharmonicity can be significantly 
higher than that of transmon qubits, allowing for faster (even subnanosecond 51 ’ 52 ) 
control pulses and reduced frequency crowding in multi-qubit systems. The addition of 
a high-quality-factor shunt capacitance to the flux qubit, together with a reduced qubit 
persistent current, has enabled us to achieve values of 7) as high as 55 ps at the qubit 
flux insensitive point. We are able to account for measured 7) values across 22 qubits 
with a single model involving ohmic charge noise, 1// flux noise, and the Purcell 
effect, with temporal variation in 7j explained by quasiparticle tunnelling. Based on this 
model, we anticipate further design optimization leading to even higher coherence will 
be possible. Finally, we used spin-locking to directly measure the photon shot noise 
spectral density, and we verified its functional fonn using a CPMG pulse sequence to 
reach a 7) of 85 ps—limited by 27)—at the flux insensitive point. These measurements 
identify photon shot noise as the dominant source of the observed dephasing, and have 
direct implications for any qubit in which the readout involves its transverse coupling to 
a resonator. 

The role of high-frequency 1// flux noise in qubit relaxation is intriguing. Our 7) data 
and their frequency dependence across 22 different qubits strongly support the 
conclusion that 1// flux noise contributes to qubit relaxation up to at least 3 GHz in 
our devices. Above 3 GHz, there is some ambiguity between ohmic flux and ohmic 
charge noise, and clarifying the roles of these respective noise sources is the subject of 
future work. A detailed understanding of such a broadband 1// -flux noise mechanism 
and its transition from classical to quantum behaviour is of great practical interest and 
awaits theoretical explanation. 
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Figure 1 | C-shunt flux qubit. 

(a) Optical micrograph of the 2.5 x 5.0 mm 2 chip, aluminum (black) on sapphire 

substrate (white, where the aluminum has been etched away), featuring two qubits (A 
and B) and a A/2 coplanar waveguide resonator ( coJIk = 8.27 GHz). The scale bar 
corresponds to 0.5 mm. (b) SEM image of the shunt capacitor (C shA = 51 fF) for qubit 

A. Each square plate of the capacitor is 200 x 200 pm 2 . The lower plate capacitively 
couples the qubit to the A/2 resonator. The scale bar corresponds to 50 pm. (c) 

Magnified view of the shadow-evaporated qubit loop and its three Josephson junctions. 
The left junction area is smaller by a factor a A = 0.42 . The scale bar corresponds to 1 

pm. (d) Device and measurement schematic. Experiments are performed in a dilution 
refrigerator at 20 mK. A global magnetic field B provides a magnetic flux bias < f ) b to 

each qubit. A qubit drive tone ( (Q d ), readout tone ( 0 ) m ), and externally generated noise ( 
P aAA , see Figs. 5 and 6) enter the A/2 resonator defined by capacitances C in andC out . 
The resonator is capacitively coupled (C g A/B ) to qubits A and B. The qubit junctions 
(“x”) have internal capacitance, C LA . B and a A . B C I A/B , and are externally shunted by 
capacitance C shA/B . Each qubit loop supports a circulating persistent current I pA/B . 
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Readout signals at the resonator output pass three isolators are amplified at 

cryogenic and room temperatures, and subsequently detected. See supplementary online 
material for more infonnation. 
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Figure 2 | Spectroscopy and T x of two capacitively shunted flux qubits. 

(a-b) Resonator and qubit spectra vs. bias current used to induce the global magnetic 
field B. The qubit transition frequencies have minima A a /2/t = 4.36 GHz and 

A b /2 n = 4.70 GHz at the qubit flux-insensitive points, which are intentionally offset in 

bias current (magnetic flux) by using different qubit-loop areas. 


(c) Energy decay functions of qubits A and B measured at their respective degeneracy 

points using the inversion-recovery pulse sequence (inset). Solid lines are exponential 
fits with decay constant7j. 
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Figure 3 | Ti variation with qubit frequency and noise modeling. 

(a) Energy relaxation time T x vs. qubit frequency for Qubit C (C shc = 9 fF , 

/ c = 275 nA, A c /2 n - 0.82 GHz ) plotted with simulated T x values for individual 

(dashed lines) and aggregate (solid line) charge, flux, and Purcell noise mechanisms. 
Absence of data around 4 GHz is related to an ancillary qubit level crossing the readout 
resonator, prohibiting qubit readout, and is not a systematic issue. Qubit C is limited by 
flux noise below about 4.5 GHz. For comparison, the functional form for ohmic flux 
noise (grey dotted line) is incompatible with the data below 3 GHz; above 3 GHz, its 
role cannot be readily distinguished from charge noise (see text). Shaded region 
indicates the range of predicted T\ in the presence of 0 - 1.0 quasiparticle. 


(b) Flux noise spectroscopy performed on Qubit B using Ramsey interferometry (red) 
and Ti p spin locking (blue) to detennine parameters = (l .4 /H> () )' /Hz and y - 0.9 
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for the inverse frequency flux noise (black dashed line) for qubit C (panel a). Green and 
black dots: inferred ohmic flux noise S 0 based on measured 7, in panel a. 

(c) Energy relaxation time T x vs. qubit frequency for Qubit B (C sh B = 51 fF , 

/ B = 49 nA, A b /2 n - 4.70 GHz). T x is sensitive predominantly to ohmic charge 
noise within 5-6.5 GHz range. Scatter in 7, is attributed to quasiparticle fluctuations. 
Cluster of lower T ] values near 5.5 GHz is due to interaction with the f X2 transition. 

Shaded region indicates the range of predicted T\ in the presence of 0 - 1.0 
quasiparticle. 

(d) T x -values for 22 qubits with widely varying design parameters, measured at their 
degeneracy points, and plotted against predicted T x values (dashed line) determined 

from numerical simulations using a single model with fixed noise levels (see main text). 
Practically indistinguishable data points (eight in total) are indicated with arrows. 
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Figure 4 | Ti temporal variation and quasiparticles. 

(a) Energy relaxation measurements (Set 1 and Set 2) at co L] = A B for qubit B. Each set 

comprises the average of 50 individual decay traces acquired sequentially in four- 
minute intervals. Set 1 exhibits purely exponential decay with JJ = 55 jus , whereas set 2 

(acquired 17 hours after set 1) exhibits a non-exponential decay function. The black line 
is a fit to Eq. (3) assuming the non-exponential decay is due to quasiparticle fluctuations 
(see text). Inset: tabulation of the values obtained from fitting functions. The * indicates 
an assumed value from set 2 (not a fit value). 

(b) Histograms of T ] values obtained by exponential fits of the individual traces forming 
the two data sets in panel (a). For set 2, the fitting is restricted to the first 40 jus to 
capture primarily the fast initial decay. 
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(c) Energy relaxation measurements at co q = A c for qubit C. The exponential decay 
function is manifest as a linear fit on the log plot with time constant T x - 0.49 jus . 

(d) Histograms of Rvalues obtained from repeated measurements of qubit C. Both the 
exponential decay function (panel c) and the consistently tight T x distribution (panel d) 

indicate a relative insensitivity to quasiparticle number fluctuations. 
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(a) 2y /2^=0.9 MHz 






Figure 5 | Calibration of engineered noise. 

(a) Resonator transmission spectra measured with the qubit prepared in the ground and 
excited states. In contrast to typical transmon qubits with co q < oo r , an excited-state C- 

shunt flux qubit shifts the resonator to higher frequencies because of interactions with 
higher-level qubit transitions. 

(b) Engineered thermal photon noise source. A coherent tone near the resonator 
frequency is mixed with white noise of nominally equal power from two independent 
arbitrary waveform generators (AWGs) applied to the in-phase (I) and quadrature (Q) 

ports of the I/Q mixer. The AWG noise bandwidth (80 MHz) is much greater than the 
cavity linewidth, creating effectively a thermal photon noise source with power P add . 


22 

































(c) Qubit spectral line shape (dots) and lorentzian fits (solid lines) for various added 
noise powers P add . The equivalent photon population n aAA added to the resonator is 

derived from (d). The blue trace corresponds to no added noise from the source in (b). 

(d) Stark-shifted qubit frequency vs. applied noise power (dots). Colored dots 
correspond to traces in (c). Combining the linear fit (solid line) with the first-order 

dependence of Stark shift on photon population (Eq. 4) yields the power-per-added- 
photon dn aAd /dP aA& = 1/(17 //W) in the low-power limit. Inset: wider range of applied 

noise powers; dashed box indicates the range in the main panel. At large photon 
populations ( n aid > 1) the frequency shift becomes nonlinear, following Eq. (43) in Ref. 

40 (solid line). 

(e) Spin-echo decay (dots) with exponential fit (solid lines) for several values of added 
photons. Inset: spin-echo pulse sequence. 

(f) Spin-echo pure-dephasing rate (echo decay rate without the T x contribution) plotted 

vs. injected photon population (dots). The linear fit (solid line) has slope 
d^pt/dn^ = 2.6xl0 6 s” 1 , in agreement with the value of 2.5xl0 6 s _I calculated from 

Eq. (5). The intercept indicates a residual photon population n 0 = 0.006 in the 
resonator. 
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Figure 6 | PSD of photon fluctuations in the resonator. 

(a) Noise power spectral densities (PSDs) extracted from spin-locking ( T lp ) relaxation 
experiments (Inset, top-left: T lp pulse sequence) measured for different locking (Rabi) 
frequencies (0.1-100 MHz) and added noise photons n aAd . Colored dashed lines 
indicate expected lorentzian noise spectra [see Eq. (6)] assuming a constant white-noise 
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offset. The vertical grey line indicates the 3-dB point of the lorentzians, coinciding with 
the resonator decay rate k/ 2n =1.5 MHz . Inset, bottom-left: T lp decay traces for 

different photon populations at fixed locking (Rabi) frequency 

<u/2;r = Q R /2;r = 500 kHz . Inset, top-right: T x decay traces at different locking (Rabi) 
frequencies for n aii =0.115 . 

(b) Decay times T 2Cm( , for the Carr-Purcell-Meiboom-Gill (CPMG) sequence measured 
vs. number of n -pulses, N n , with n add = 0 (no added noise photons). The CPMG pulse 

sequence (inset) acts as a band-pass noise filter centred at a frequency proportional to 
N k through the pulse spacing Ar (see main text). At N, = 200 , the filter frequency 

approximately equals the cavity decay rate (dashed line). For N n < 200 , the filter 

samples the flat low-frequency portion of the lorentzian PSD, yielding a constant decay 
time r 2CPMG ~ 40 //s . For N„ > 200, the filter traverses the roll-off region of the 

lorentzian. As the sampled noise decreases, the decay times increase, approaching the 
limit set by energy relaxation (f 2CPMG « 27^) for A t > 1000. 
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Supplementary Figure 1: Process flow schematics of key steps for capacitively-shunted flux qubit fabrication, a. 

Preparation of MBE aluminum (red) on outgassed C-plane sapphire substrates (gray), b, Patterning of the MBE aluminum 
into the shunt capacitor (representative square shunt capacitor geometry is shown), resonator center line, and surrounding 
ground plane, c, Patterning of the aluminum qubit loop (yellow), which contains three aluminum Josephson junctions. The 
loop contacts the shunt capacitor as illustrated. 
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Supplementary Figure 2: Experimental measurement schematic, a. Diagram of the input-signal generation of readout (top 
panel, cyan), qubit drive (mid panel, magenta) and noise (bottom panel, brown) at room temperature. The signals are combined 
before being sent to the fridge, b. Relative timing of the signals generated in a. The dashed lines indicate the window within 
which the readout pulse is turned on while the noise is turned off. The representative qubit drive is a spin-echo pulse sequence, 
two 7r/2-pulses and a midpoint 7r-pulse. The actual pulses assume a Gaussian envelope, c. Wiring inside the refrigerator. 
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Supplementary Figure 3: Circuit diagram of the C-shunt flux qubit. Josephson junction 3 is a times the size of junctions 1 
and 2. The external flux / e is defined as the magnetic flux threading the loop formed by the three junctions. Nodes A and B 
represent the superconducting islands. The shunt capacitor has the dominant capacitance in the circuit. Node A is also 
capacitively coupled to a superconducting resonator (not drawn). The green dashed box highlights the part of the circuit 
resembling the conventional persistent-current flux qubit. 
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Supplementary Figure 4: Potential profile of the C-shunt flux qubit, a. A contour plot of the potential energy U{p\ , <^ 2 ) at 
f b = 0 for a = 0.4 (C-shunt flux qubit, on the left) and a = 0.8 (traditional flux qubit, on the right). In the C-shunt case, the 
square-shaped area centered around (p i,pi) = (fci -27r, /c2-27r) marks the single well. In the traditional case, the 
figure-eight-shaped area marks the double well. The gray arrow indicates the <^ m = p>\ — pi direction, b. The potential 
function along the p m direction at pi + pi = 0 at f b = 0 and f b = 0.03 for both a = 0.4 and a = 0.8 cases. 
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Supplementary Figure 5: Circuit schematic for the C-shunt flux qubit and its reduction for use in simulation, a. 

Schematic of the qubit-resonator system and the islands / nodes used in the simulation. For clarity the labeling follows that 
used for cQED transmons in Ref. 1. b. Equivalent circuit schematic with islands 1 and 5 grounded. Filled black circles 
indicate the five nodes that will be used in the full simulation, c. First reduction of the circuit in (b) containing four nodes 
(labeled in blue). See text for details, d. Second reduction of the circuit in (b) containing three nodes (labeled in blue). An 
effective shunt capacitance C^h accounts for the capacitances Ct, Cb and C a . See text for details. 
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Supplementary Figure 6: High-power spectra of qubits A and B with simulations to match the higher-level transitions. 

The horizontal axis corresponds to the current in the “bobbin” coil of wire local to the qubits to apply a magnetic flux bias. 
Solid lines are simulation results that give the best match to the visible 0-1, 0—2 (two photon), 1—2, 0—3 (two photon) and 0—2 
transitions. A precise matching to the experimental spectra required allowing these parameters to vary somewhat from their 
design values (Supplementary Table 1). Nonetheless, the simulation parameters (e.g., J c , C a h, etc.) are generally same for the 
two qubits, with the primary exception being the junction sizes. 
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Supplementary Figure 7: Two-tone spectroscopy of qubit B measured at the optimal bias point ('[>,. = <f> 0 /2). The data 
are taken by first applying a 20 ns microwave pulse at = cjoi/27t = 4.701 GHz to drive the qubit to the |1) state, 
followed by a long (50 fis) low-power microwave tone. The frequency of the second tone is swept to perform spectroscopy. 
The resulting spectrum shows two peaks, one at 4.701 GHz that corresponds to driving the 1—0 transition, and one at 
5.191 GHz that corresponds to the 1—2 transition. The measured anharmonicity is (cui 2 — o>oi)/27r = 490MHz. 
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Supplementary Figure 8: Circuit schematic used to account for electric and magnetic noise coupled to the qubit. See 

main text Supplementary Note Supplementary Note 8 for details. 
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(a) Free- and driven-evolution analogy 


Free evolution Driven evolution 

Qubit frame Rotating frame 




*t 


(c) Dynamics in SL-3 Seq. 



Supplementary Figure 9: Spin locking and its implementation for noise spectroscopy, a. Analogy between free- (left) and 
driven-evolution (right) dynamics. The free evolution is described in the qubit frame {x, y, z} while the driven evolution is 
described in the rotating frame {X, Y, Z}. The two cases differ in the orientation and size of the static field, i.e., “tVqCTz/2" 
(red arrow) versus “wrux^" (yellow arrow). The corresponding quantization axes (green arrows) and 
longitudinal/transverse relaxations (gray arrows) are defined with respect to the static field, b. Standard three-pulse 
spin-locking sequence (SL-3). The long driving pulse is 90°-phase-shifted from the 7r/2-pulses, and its length r is the scanned 
parameter to record the rotating-frame relaxation, c. Bloch sphere representation of the rotating-frame qubit dynamics under 
SL-3. The purple arrows represent the polarization of the qubit states, while the magenta arrows indicate the driving-field 
orientation. The qubit is initially prepared in its ground state (I). The first 7r/2-pulse rotates the qubit by 90° into the 
equatorial plane (II). The second 90°-phase-shifted continuous driving pulse, of duration r, is then aligned with the qubit 
state, effectively locking the qubit along X. During the pulse, the qubit undergoes relaxation in this rotating frame towards its 
steady state (III). The final 7r/2-pulse projects the remaining polarization onto Z (=z) for readout (IV). 
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Left column figures: 

Model includes 
ohmic charge noise 
(same as main text) 


Right column figures: 

Model replaces 
ohmic charge noise 
with ohmic flux noise 




Qubit frequency (GHz) 


Qubit frequency (GHz) 





Supplementary Figure 10: Comparison of a model using ohmic charge noise with one using ohmic flux noise. Panels a. c, 
and e reproduce the panels in Fig. 3 from the main text, and the simulation traces contained therein are derived from a model 
involving ohmic charge noise, 1// flux noise and the Purcell effect. For comparison, panels b, d, and f replace the ohmic 
charge noise with ohmic flux noise. The agreement with the data is similar between the two models, making it challenging to 
distinguish between ohmic charge noise and ohmic flux noise. It is only in panels a and b (device C) at frequencies above 6 
GHz where the qubit is especially sensitive to ohmic charge noise that a plausible distinction can be observed. 
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Supplementary Figure 11: System filter functions corresponding to representative control sequences. Filter functions of 
Ramsey, spin-echo and CPMG (N = 2,4,8,16) sequences, assuming r = 1 fis in all cases. 
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Sample Shunt-cap. and Coupling Junction Parameters Qubit Parameters at <3>b = 0 Measured 



main 

cap 

gap 

C Bh 

/3c 

a 

Jc 

Eja 

Eca 

A 

A 

Veg 

h 

Tf 

t q 

Tf 

jisim 

'pmeas 

A 

# 

text 

type 

(gm) 

(fF) 



(fiAlfim 2 ) 

(GHz) 

(GHz) 

(GHz) 

(GHz) 

(mV) 

(nA) 

(Ms) 

(M s ) 

(M s ) 

(MS) 

(M s ) 

(GHz) 

1 

IDC 

40 

21.3 

0.19 

0.55 

6.7 

139 

0.74 

3.36 

6.15 

4.28 

194 

4.4 

21 

742 

3.6 

4.0 

4.4 

2 

IDC 

40 

21.3 

0.19 

0.55 

6.7 

139 

0.74 

3.36 

6.15 

4.28 

194 

4.8 

19 

742 

3.8 

2.6 

4.9 

3 

IDC 

20 

21.3 

0.19 

0.59 

6.6 

127 

0.75 

1.50 

8.37 

2.46 

206 

1.8 

155 

3014 

1.8 

1.6 

1.8 

4 

IDC 

20 

17.5 

0.19 

0.63 

6.6 

137 

0.88 

0.88 

11.65 

1.66 

275 

0.5 

821 

2894 

0.5 

0.5 

0.7 

5 

c 

IDC 

40 

9.0 

0.31 

0.68 

6.6 

133 

1.46 

0.82 

17.8 

1.76 

275 

0.5 

685 

2894 

0.5 

0.5 

0.8 

6 

IDC 

40 

21.3 

0.19 

0.55 

6.6 

137 

0.74 

3.36 

6.10 

4.30 

191 

2.7 

37 

742 

2.5 

1.9 

2.5 

7 

IDC 

40 

22.7 

0.31 

0.53 

5.9 

93.1 

0.73 

3.88 

4.40 

4.74 

126 

8.4 

22 

241 

5.9 

4.0 

3.5 

8 

IDC 

20 

27.6 

0.21 

0.52 

5.9 

87.3 

0.62 

3.94 

3.51 

4.45 

113 

15 

16 

408 

7.6 

7.0 

5.4 

9 

IDC 

20 

17.9 

0.23 

0.56 

5.9 

105 

0.88 

3.31 

6.29 

4.61 

152 

6.3 

21 

536 

4.8 

5.0 

3.9 

10 

IDC 

20 

17.9 

0.23 

0.55 

5.9 

105 

0.88 

4.19 

5.53 

5.34 

145 

6.7 

16 

305 

4.7 

5.0 

3.7 

11 

IDC 

40 

21.3 

0.19 

0.55 

5.9 

122 

0.74 

3.35 

5.75 

4.30 

171 

4.4 

28 

762 

3.7 

2.3 

3.3 

12 

IDC 

40 

21.3 

0.19 

0.54 

5.9 

122 

0.74 

4.14 

5.05 

4.90 

163 

5.8 

18 

467 

4.3 

3.5 

4.1 

13 

IDC 

20 

26.8 

0.15 

0.52 

3.0 

44.3 

0.64 

3.19 

2.73 

4.05 

61.8 

32 

33 

1378 

16 

12 

3.2 

14 

IDC-R 

10 

51.4 

0.14 

0.43 

2.4 

36.2 

0.35 

3.94 

0.88 

3.43 

43.9 

82 

34 

643 

23 

23 

4.2 

15 

IDC-R 

10 

51.4 

0.14 

0.43 

2.4 

36.2 

0.35 

3.94 

0.88 

3.43 

43.9 

82 

34 

643 

23 

22 

4.2 

16 

IDC-R 

30 

51.9 

0.14 

0.43 

2.4 

36.2 

0.35 

3.94 

0.87 

3.43 

43.9 

82 

34 

643 

23 

24 

4.2 

17 

IDC-R 

30 

51.9 

0.14 

0.43 

2.4 

36.2 

0.35 

3.94 

0.87 

3.43 

43.9 

82 

34 

643 

23 

22 

4.2 

18 

SQR 

40 

51.0 

0.089 

0.43 

2.9 

43.0 

0.35 

4.33 

0.88 

3.60 

50.3 

68 

147 

408 

42 

35 

4.7 

19 

SQR 

40 

51.0 

0.089 

0.43 

2.9 

43.0 

0.35 

4.68 

0.88 

3.74 

50.1 

75 

124 

330 

41 

44 

5.1 

20 

SQR 

40 

51.0 

0.089 

0.42 

2.4 

36.2 

0.35 

3.96 

0.80 

3.45 

43.9 

79 

184 

567 

50 

42 

4.0 

21 

A 

SQR 

40 

51.0 

0.089 

0.43 

2.4 

36.2 

0.35 

3.96 

0.91 

3.45 

43.9 

85 

169 

567 

52 

44 

4.4 

22 

B 

SQR 

40 

51.0 

0.089 

0.42 

2.4 

36.2 

0.35 

4.29 

0.83 

3.60 

43.5 

92 

145 

474 

50 

55 

4.7 


Supplementary Table 1 : Parameters for 22 qubits studied in this work, spanning shunt capacitance values C s h = 9... 52 fF. Color-highlighting corresponds to the 
samples A, B and C presented in detail in the main text. All parameters are design values or are derived from simualtion, except for the critical current density J c and 
the measured relaxation time T™ eas . Three capacitor types were studied: IDC is “interdigital capacitor”; IDC-R is an IDC with slightly rounded comers (negligible 
when compared with actualized standard IDCs, but included here for completeness); and SQR is square capacitor shape. Gap is the spacing between capacitor 
features ( e.g IDC fingers). 0 C quantifies the coupling strength between the qubit and the resonator. Ej a and Eca are respectively the Josephson energy and charging 
energy of the small junction. A is the qubit frequency at the qubit flux-insensitive (degeneracy) point <3?b = 0. A is the anharmonicity. V eg is the simulated island 
node voltage (related to the charge transverse matrix element). I p is the simulated persistent current (related to the flux transverse matrix element). Tp, T®, and T\ 
are respectively the simulated Ti values due to flux noise, charge noise, and the Purcell effect (all at <3?b = 0). Tf im is the net simulated T\ value due to these 
processes, and and T™ eas is the measured relaxation time at <3>b = 0. 

































Supplementary Note 1 Materials and Fabrication of the High-Q Capaci¬ 
tor 

The capacitively-shunted flux qubits studied in this work were prepared using the following steps: 

1. Growth and patterning of high-quality-factor (high-Q) aluminum films using molecular beam epitaxy (MBE). 

2. Patterning and evaporation of the superconducting qubit loop and Josephson junctions. 

3. Dicing and packaging. 


Supplementary Note 1.1 Growth and patterning of high-Q aluminum 

High-Q aluminum films were deposited on 50-mm C-plane sapphire wafers in a Veeco GEN200 MBE system 
with a growth chamber base pressure of 10 _n torr. The wafers were cleaned in piranha solution (sulfuric acid and 
hydrogen peroxide) prior to loading into the MBE system. The wafers were annealed in the MBE system at 900 °C 
to facilitate outgassing and sapphire surface reconstruction, after which 250 nm of aluminum was deposited at a 
growth rate of 0.025 nrn/s and a substrate temperature of 150 °C (Supplementary Figure la). 

The high-Q aluminum was patterned using contact lithography and wet-etched using Aluminum Etchant - 
Type A (Transene Company, Inc.) into the following device features (Supplementary Figure lb): shunt capacitors, 
coplanar waveguide (CPW) resonators, ground planes, and optical alignment marks. 

Supplementary Note 1.2 Patterning the qubit loop and Josephson junctions 

The qubit loop and junctions were formed using double-angle evaporation of aluminum through Dolan-style 
bridges (Supplementary Figure lc) [2]. The free-standing bridges were realized using a bilayer mask compris¬ 
ing a germanium hard mask on top of a sacrificial MMA/MAA layer [MicroChem methyl methacrylate (MMA 
(8.5)/MAA EL9)]. The qubit loop and junctions were patterned using electron-beam lithography (Vistec EBPG5200) 
and ZEP520A resist (ZEONREX Electronic Chemicals). This pattern was transferred into the Ge layer using a CF 4 
plasma, and the underlying MMA/MAA resist was under-etched using an oxygen plasma to create free-standing 
bridges. Prior to the aluminum evaporation, an in situ argon ion milling was used to clean exposed contact points on 
the MBE aluminum to ensure superconducting contact with the evaporated aluminum. The qubit loops and junc¬ 
tions were realized with two separate angle-evaporated aluminum layers; between the two aluminum evaporation 
steps, static oxidation conditions were used to prepare junctions with a certain critical current density. 

Supplementary Note 1.3 Dicing and packaging 

Devices were diced into 2.5 x 5 mm 2 chips (as shown in the manuscript in Fig. la) that were mounted into gold- 
plated copper packages. Aluminum wirebonds were used for both signal and ground connections between the 
device and package, as well as to connect the ground planes of the CPW resonator to prevent slotline modes. 
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Supplementary Note 2 Measurement Set-up and Protocol 

We performed our experiments at MIT in a Leiden Cryogenics (CF-450) dilution refrigerator with a base temper¬ 
ature of 15 rnK. The device was magnetically shielded with a superconducting can surrounded by a Cryoperm-10 
cylinder. All electrical leads were attenuated and/or filtered to reduce noise. 

Supplementary Note 2.1 Outside the dilution refrigerator 

The electronic setup for generating the readout-pulses, control-pulses, and artificial photon noise is shown in Sup¬ 
plementary Figure 2a. All time-sensitive instruments are synchronized with a Stanford Research Systems FS725 
Rubidium Frequency Standard. The readout-pulse is generated by an Agilent 8267D Vector Signal Generator 
(PSG), gated by a Tektronix 5014b arbitrary waveform generator (TEK). This pulse is typically a few microseconds 
long, and the tone is set at the resonator frequency dressed by the qubit in its ground state (cuf/27r). The control- 
pulse envelope for driving the qubit is generated by the TEK and mixed with a qubit-frequency tone (u> q /27r) from 
a second PSG using its internal I-Q mixer. These pulses are further gated to reduce unwanted mixer leakage. 
Artificial photon noise is generated by up-converting (quadrature-mixing) 80-megahertz-wide white-noise signals 
from two Agilent 33250A arbitrary waveform generators (AGI) with a tone near the cavity frequency from a third 
PSG. We confirmed that the noise source behaved consistently using either one port or both ports of the I-Q mixer. 
For example, with equal powers applied to each port, the output is simply a doubling of the photon noise power 
generated from a single port. The carrier frequency for the noise is chosen to be ~10 MHz away from the readout 
frequency, much larger than the cavity linewidth. In addition, we gate the noise off during the readout pulse using 
the TEK. The three signals are combined (relative timing illustrated in Supplementary Figure 2b) and sent to the 
input port on the refrigerator. This sequence constitutes a single experimental trial, and it is repeated typically 
10,000 times with a period of 200-400 /rs. 

Supplementary Note 2.2 Inside the dilution refrigerator 

The device is enclosed in a copper package, which is itself mounted inside an aluminum box to shield the device 
from external electromagnetic radiation and magnetic field fluctuations (Supplementary Figure 2c). A small coil 
antenna is mounted under the package lid and is used to bias the qubit with a static magnetic field (not shown). A 
Yokogawa 7651 dc source provides the bias current to the coil using twisted pair wires passing through an RC low- 
pass filter with cutoff around 100 kHz. On the input side, there is a total attenuation of 49 dB arising from discrete 
attenuators (XMA Corporation) at various temperature stages. In addition, there is a relatively small amount of 
distributed attenuation due to loss in the coaxial cables. On both the input and output sides, a high-pass filter (RLC 
F-18948, 4 GHz cutoff) and a low-pass filter (RLC L-3615, 12.4 GHz cutoff) provide a net 4-12.4 GHz passband. 
After the output filters, 3 isolators (Quinstar / Pamtech, model CWJ1019KS414, 3-12 GHz, with approximately 
15-20 dB isolation each) are mounted on the mixing chamber. The output signal is amplified by a JPL/Caltech 
cryogenic preamplifier (1-12 GHz, 30 dB gain). The output port outside the refrigerator is followed by a room- 
temperature amplifier (MITEQ, AMF-5D-00101200-23-10P, 0.1-12 GHz, 43 dB gain) before mixing with a local 
oscillator (LO) tone (w ro /27r — wlo/27t = 50 MHz ) for heterodyne detection (not shown). After the mixers, 
the signal is digitized using an Acquiris U1084A analog-to-digital converter (ADC) and digitally demodulated to 
extract the amplitude and phase of the readout signal for qubit-state estimation. 
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Supplementary Note 3 Two-Level-System Model of C-Shunt Flux Qubit 

Supplementary Note 3.1 Parameterization of the two-level-system Hamiltonian 

In the main text, we elect to parameterize the C-shunt flux qubit using the familiar notation of the conventional 
flux qubit [3], albeit with important generalizations due to the influence of higher energy levels. In the laboratory- 
frame, the two-level-system Hamiltonian for an individual C-shunt flux qubit near flux-degeneracy and coupled to 
a CPW resonator is: 

TL = ^ [A ($ b ) <5x + £ ($b) o z ] + huj r (a^a + 1 /2) + hg (<h b ) <5' y (d t + a). (1) 

Here, the three terms are respectively the qubit, resonator, and qubit-resonator interaction Hamiltonians, <r x yjZ are 
the Pauli operators for the qubit, <l>b = 3> e — ^o/2 is the flux bias due to an external magnetic flux <l> e (T>o = h/2e 
is the superconducting flux quantum), a^(a) is the raising(lowering) operator for resonator photons, oj t is the 
resonator angular frequency and g is the qubit-resonator coupling angular frequency. In this coordinate system, the 
qubit shunt-capacitor couples transversally to the qubit through a y . Similarly, charge fluctuations connect through 

(Ty. 

Within our parametrization, the circulating current states of the C-shunt flux qubit have energies e(fl>b)/2 = 
±7 n i(fl’b) fl’b? where 7 m (<£>b) = [7 p i (^b) — P p o (^b)]/2 is related to the difference in the flux-dependent circulat¬ 
ing (persistent) currents /pc/fl’b) and 7 p i(<l>b), and these states hybridize with a flux-dependent energy fi.A(fl>b). 
With this generalization, a flux-dependent offset energy he 0 s (<3>b) = lots (^b) fl’b is removed, where 7 0 ff ($b) = 
[Ipi (*b) + Ipo ($b)]/2 ■ In addition to their flux dependence, the currents 7 p o, p i(<l , b) generally do not have the 
same magnitude, nor need they have opposing sign. This is in marked contrast to the conventional flux qubit, 
where 7 p o = —7 p i = 7 p with both 7 p and A being independent of flux within a certain region about $b = 0. The 
fact that 7 p o,pi (d> b ) and A ($b) are flux-dependent for the C-shunt flux qubit reflects the non-negligible role of its 
higher energy levels. Consequently, we numerically diagonalize the full qubit-resonator Hamiltonian for parameter 
extraction and modeling. 

The extent to which higher levels of the qubit influence the parameters in Eq. (1) is determined by the ratio 
/cu q where u> q = -y/ £ («I>b ) 2 + A($b ) 2 is the qubit frequency, and uj p fc) are the plasma frequencies of the 
additional oscillator-like modes of the circuit. For conventional flux qubits, these modes are typically much higher 
in energy than the qubit excited state and can be neglected. However, as the shunt capacitance is increased, at least 
one of these modes shifts to lower frequency and can eventually become comparable to the qubit frequency itself. 
For example, one of these modes is the small junction “plasma frequency” w PjQ = 1/ \jLj, a C a determined by 
its Josephson inductance Lj Q and its total capacitance C a (including the shunt capacitor). This is the frequency 
of the oscillator-like degree of freedom for the qubit in each of the two wells of its double-well potential profile 
(each well is associated with a circulating persistent current 7 p o, p i), and this frequency is ideally far above u/ q . 
For a conventional flux qubit, w p /w q > 10, and the parameters in Eq. (1) are largely flux-independent over 
a wide range of flux bias about fl>o/2 (see Ref. 4 for experiments that study the flux range required to access 
higher levels for a particular conventional flux qubit). As the ratio w p /co q is decreased, the excited states of 
the plasma mode shift closer to the two qubit levels and quantum mixing (hybridization) occurs, resulting in 
increased quantum fluctuations of the current. Correspondingly, there is a decrease in the range in flux over 
which the parameters in Eq. (1) are essentially flux independent. Eventually, when /u> q > 1, the qubit can be 
descried as a weakly-harmonic oscillator (much like a transition), where the only remnants of the persistent currents 
are small, state-dependent displacements of the oscillator current from zero, which are much smaller than their 
quantum fluctuations. Typically, the resulting weak anharmonicity is inverted relative to that of the transmon. The 
highest-coherence C-shunt flux qubits considered in this work are in an intermediate regime, having o> p fc ' /oj q > 2. 
Their anharmonicity is around 500-900 MHz for the longest lived devices, generally larger than the 200-300 MHz 
observed in transmons. 

The role of higher energy levels is described in Section Supplementary Note 4, along with a comparison 
between the conventional and C-shunt flux qubits. Ultimately, to account for the higher energy levels, the two- 
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level parametrization must be deduced from experiment, an analytic treatment (Section Supplementary Note 4), or 
simulations (Section Supplementary Note 5). 


Supplementary Note 3.2 Dispersive two-level-system Hamiltonian 

The unitary U\ = exp [i (| — #)<x y /2] rotates Eq. (1) from the laboratory frame to the qubit frame, 

U\HU\ = Hojqfrz/2 + huJria)a + 1 / 2 ) + hga y (a * + a) , ( 2 ) 

where 6 = arctan(e/A). The transverse qubit-resonator coupling is unaffected by the transformation, since the 
flux bias 4>b rotates the qubit quantization axis within the x—z plane (laboratory frame) and the capacitive coupling 
is transverse to this plane. 

In the dispersive regime, |^ q r | = |w q —w r | g, after a second unitary transformation U 2 = exp [—+ 
a t <T —)], the Hamiltonian can be approximated: 

'Wdisp = U 2 U 1 HU 1 U 2 ~ Sw q «r z /2 + hjj T (a)a + 1/2) + H2x(a^a + l/2)o z /2 . (3) 

Here, w q = V 7 ^(^b ) 2 + A ($ b ) 2 is the qubit angular frequency and x(fl’b) is related to the qubit-state-dependent 
dispersive shift 2x($b) = g 2 (<I > b)/<5q,r( < I > b) of the resonator used for readout. The last term includes the Stark 
shift Agtark = 2x(*£*b) n due to the resonator photon number h = aid and the Lamb shift ALamb = x(^b) due 
to the resonator zero-point energy. Higher levels of the qubit generally play an important role in constructing the 
flux-dependent dispersive shift. While it is straightforward to measure this experimentally, one must go beyond a 
two-level approximation to calculate it accurately. 


Supplementary Note 4 Analytic Treatment and Comparison of the Con¬ 
ventional and C-Shunt Flux Qubits 

Supplementary Note 4.1 Conventional persistent-current flux qubit 

Before developing the system model for the C-shunt flux qubit, we review the conventional persistent-current flux 
qubit as presented in Ref. [3]. The conventional flux qubit circuit is illustrated in Supplementary Figure 3 (dashed 
box). The qubit loop is interrupted by three Josephson junctions, and ipi ( i = 1, 2, 3) are the associated gauge- 
invariant phase differences. Two of the junctions have the same critical current I c and junction capacitance C, 
and thus the same Josephson energy Ej = 7 c fl>o/27r and charging energy Ec = e 2 /2 C. The third junction is 
smaller in area by a factor a, resulting in a reduced critical current al c and reduced junction capacitance aC. The 
corresponding energy scales are aEj and Ec/a. 

The flux quantization condition gives </?i — g>i + <£>3 = 2nf e , where f e = < T e /^ , o is the external magnetic 
flux threading the loop in units of the superconducting flux quantum $0 = h/2e. When a > 0.5 and / e « 0.5, 
the potential energy of the qubit assumes a double-well profile. The wells are associated with clockwise and 
counterclockwise circulating currents tunable by the applied magnetic flux. These diabatic circulating-current 
states tunnel-couple with a strength depending on the height of the inter-well barrier or the scale factor a. 

Tuning / e tilts the double well potential. In the vicinity of / e « 0.5, the circulating currents are of opposite sign 
and essentially equal magnitude. Consequently, in this limit, higher energy levels play little role and the two-level 
approximation is a good approximation to the full Hamiltonian. However, even for the conventional flux qubit, 
higher energy levels become important for flux biases far from f e « 0.5, where the double-well potential is tilted 
to such a degree that its higher levels influence the circulating currents [4,5], The region about f e ss 0.5 for which 
the two-level system is a good approximation is reduced as a decreases. 
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Supplementary Note 4.2 C-shunt flux qubit: quasistatics 

The C-shunt flux qubit has three features that are distinctly different from the conventional design: (i) lower 7 C , 
(ii) lower a, typically a < 0.5 (defined for three-junction qubits) and (iii) an additional large capacitor C s h = (C 
(( 1) shunting the smaller junction. As we show below, both (i) and (ii) lead to reduced sensitivity to flux noise, 

while (iii) reduces sensitivity to charge noise. 

Following a recipe similar to that presented in Ref. [3], the three-junction capacitively-shunted flux qubit can 
be described by a Hamiltonian consisting of the kinetic and potential part: 


n = f+u, 

T=\( Q + qfC-'fQ + q) , 

U = Ej{2 + a — cos(/3i — cos<p 2 — acos(2nf e + ip i — p> 2 )} ■ 


In the kinetic part, Q and q are the charges and induced charges on the islands: 


Q = —i 2e 


d 

dip i 
d 

dip2 


and q = 



( C + l + a -(C + «) \ 
\ —(C + a ) C + l + a ) 


(4) 


(5) 


The kinetic energy represents the total electrostatic energy stored in the capacitors, and is dominated by the shunt 
capacitor, since ( 1, a. The shunt capacitor largely reduces the effective charging energy, causing the system 

less sensitive to charge fluctuations. Compared with the circuit model in Ref. [3], there is no explicit gate electrode 
in our circuit. The induced charge, however, is included for modelling the charge noise, and will be discussed later 
in this section. In the static case, q A = q B = 0. The potential energy sums up the Josephson energy stored in the 
junctions. It assumes a two-dimensional periodic profile. When a < 0.5, as in this C-shunt design, there is only 
one well in each unit cell (see Supplementary Figures 4a and 4b), a distinction from the double-well profile in the 
conventional case. This leads to a smaller circulating current for the C-shunt flux qubit, reducing its sensitivity to 
flux fluctuations. 

By choosing < p p = (ip i + p>2)/‘Z and <Pm = (‘Pi — <F2)/2 as coordinates and using the Cooper-pair number 
operators h a = —i d/dip a ( cr = p, m), we have the reduced Hamiltonian 


T-L = -E c , P h 2 p + -B c ,mnm + e j{ 2 + a - 2 cos ip p cos ip m - a cos(2irf e + 2<p m )} , 


( 6 ) 


where Eq, p = e 2 /{C /2 ) and Ec. m = e 2 / C(( + a + 1/2) are the effective charging energy for the p-mode and m- 
mode respectively. Ideally, because the introduction of the shunt capacitor ensures that Eq v « 2(Ec,m Ec, m , 
we can safely omit the p-mode, since the characteristic frequency H p = ^/Ec, p Ej p /h « 2ir x 50 GHz is much 
higher. Here, Ej p = 2Ej is the effective p-mode Josephson energy. In a realistic design, one must also assess the 
role of additional plasma modes, for example, due to capacitance to ground, and their characteristic frequencies. 
The simplified Hamiltonian becomes 


— 2 


+ Ej{-2 cos ip m + a cos(27r/ b + 2p m )} 


(7) 


where fb = fe~ 0.5 is the reduced external flux bias away from one-half flux quantum. 

To understand the energy structure of the system, we Taylor-expand the potential function U = Ej{—2 cos p m + 
acos(27r/b + 2<p m )} around the well minimum 


U(p m ) = Y, U(k) ^m-<P*J k , 

fc =0 


( 8 ) 
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where U ( fc ) = (l/k\)(d k U/dp^ a \ iPai=ip - n ). For a given / b , ip^ is single-valued within the range p m e [— 7r, 7r] and 
obtainable by solving the equation U'(jp'^ B ) = 0- By definition, U (°) is a constant and U < - 1 ' 1 = 0. Therefore, we 


oo 

+ Y J U < - k) <p*(b + &) k . ( 9 ) 

k =3 

On the second line, Ej m = 2£/( 2 ' is the effective Josephson energy, ip m i = ip m > = ip m — tp^ and n m / = 
—i d/dip m '. We group the quadratic potential term with the kinetic part so that their combination Ho describes 
a one-dimensional harmonic oscillator. The last line expresses the Hamiltonian in terms of the raising (in) and 
lowering ( b ) operators of the harmonic oscillator. The oscillator frequency is 121!? = The quan¬ 

tum ground-state phase uncertainty is p z = (i?c,m/4i?j.m) 1 ^ 4 [6]. After the replacement ip m i = p z (b + w), 
the perturbation matrix V can be conveniently computed in the space spanned by the Fock states of Hi,. Note 
that U ( fc ) and <p z are all /^dependent. The eigenstates can be solved for a given order of expansion. In 

general, higher-power expansion terms introduce anharmonicity as well as modulate the harmonic frequency. For 
simplicity and without loss of generality in the approach, we restrict ourselves to the subsystem spanned by the 
first-three levels {|g), |e) , |f)}. Within a first-order perturbation analysis, we have the approximate Hamiltonian 

H = h(siPb+^PPtt>\ , ( 10 ) 


can rewrite Eq. (7) as: 


H m = Ho + V 


— ( ,y /'C.ndfin' + 


= hn£\bn+-) 


where (2 and A are respectively the /^dependent harmonicity and anharmonicity. As an important example, they 
can be derived at / b = 0, where ip^ = 0. The expansion terms then become 


(7 (2) /£j = 1 - 2 a , 

U^/Ej = 0 , 
U^/Eo = 8 °^ 
U^/Ej = 0 , 


( 11 ) 


All odd-power terms disappear due to symmetry of the potential at / b = 0. The term t/* 2 ) provides a quadratic 
potential and thus harmonicity when a < 0.5, which is consistent with the "slightly anharmonic oscillator" model 
developed here. When a > 0.5, as in the conventional flux qubit, the potential becomes a concave function with a 
rising barrier which splits the landscape into double wells. The quartic potential from is the leading term that 
introduces anharmonicity. Keeping terms up to [/A 4 ), we have hfl = \/Ec,mE]^ m + Ai?c,m and HA = \Ec, m , 
where Ej m = 2(1 — 2a)Ej and A = 8 (i^a) • Our present C-shunt designs are within the regime of a > 0.125, 
giving a positive anharmonicity in contrast to the negative one which is characteristic in transmon qubits. 
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Supplementary Note 4.3 C-shunt flux qubit: time-dependent fluctuations 

We now consider flux noise and charge noise within this analytical model. Small flux fluctuations (Sf) affect this 
anharmonic system by modulating the potential energy. That is 

SHf = 5U = —2naEjsm(2nfb + 2ip m )Sf 

= —2TraEj{shi(j) cos(2<p m /) + cos <f> sin(2<^ m / )}Sf 

« —2naEj{sm<fi (1 — 2ip 2 (b + b l) 2 ) + cosf>2ip z (b + b^)}Sf , ( 12 ) 

where <j> = <p(fb) = 2nfb + 2<Pm(/b)- The (6 + E) 2 term contains longitudinal modulation while the (b + E) 
term provides only transverse modulation, that is, a coupling matrix element between adjacent levels. 

In the spectral domain, fluctuations near frequency ft are of particular interest because of their ability to induce 
g —e transitions, and thus play an essential role in qubit driving and relaxation. Consider a coherent drive Sf oc 
cos(wdi) with driving frequency wa ~ ft. In the interaction frame, one can show that both the fast longitudinal 
oscillations and transverse counter-rotating terms are negligible in the weak driving limit. This simplifies the 
three-level interaction-frame Hamiltonian to 


^ / 2Aw + A \/2H r /2 0 \ 

H = h\ \J2 ftn/2 Aw Hr/2 

V o Hr/2 0 ) 


(13) 


where Aw = ft — Wd is the detuning of the drive from the g—e transition frequency. In the regime of weak drive, 
namely ui q ,A 3 > Hr, the system is effectively protected from transitions to the third level. Therefore, we can 
further reduce the system to two levels, 

— \ (^q + Tm$o Sf <T X ) , (14) 


where w q = w q (/b) = ft as in Eq. (10) and I m = / m (/b) = —8naip z cos</i?j/4>o. Here, I m is the //-dependent 
current difference between the parameterized circulating-current states (See section Supplementary Note 3). The 
scaling I m oc aEj^E^ 4 ^ 1 / 4 indicates the efficiency of reducing flux-noise sensitivity by lowering a and Ej, 
as we implemented in the C-shunt design. For example, qubit sample B in our device (see main text) has a « 0.4, 
Ej/h « 65 GHz and ip z « 0.28, leading to I m « 55 nA. 

On the other hand, charge fluctuations (<5</ A and Sq B ) invoke perturbation via the kinetic energy T. Expanding 
T in Eq. (4) and ignoring the p-mode (as before), we have the perturbation Hamiltonian, 


6U„ = ST = <Sq T C _1 Q 


■ [-i (S - *')] ( *9*-*>•)■ 


(15) 


On the second line, we used the transformation n = — i n z (b — w), where n z = (Tlj. m /4i?c,m ) 1 ^ 4 is the quantum 
ground-state uncertainty in Cooper-pair number. Similar to the flux-noise argument above, charge fluctuations also 
connect to the system transversely, and moreover, are also orthogonal to the flux. In addition, the perturbation 
depends only on the differential mode of the induced charges between islands A and B, the branch charge across 
the small junction. 

Adding STl q to the two-level approximated Hamiltonian in Eq. (14), we have 

if = \ < 7 Z + / m <f>o Sf < 7 X + n z E c ,mSn m < 7 y ) , (16) 


where Sn m = ( Sq A — Sq B )/(—e) is the differential electron number fluctuation. The charge-noise sensitivity is 
n z Ec,m oc E'Q A m oc C -3 / 4 . Therefore, the introduction of a large shunt capacitor makes the system less sensitive 
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to charge fluctuations. In qubit sample B (see main text), n z Ec m /h « 1.2GHz. Assuming Johnson-Nyquist 
noise from a 377-fi resistor (free-space impedance) and a parameterized gate capacitance of 0.03 C = 0.12 fF, we 
find a charge-noise-limited T\ of about 60 fis. 

The capacitive coupling of the qubit loop to the resonator can be modelled in a similar way as the node charges. 
The voltage fluctuations on the coupling capacitor C g induce charge fluctuations, and hence enter the Hamiltonian 
through the same channel as charge. Consequently, the interaction Hamiltonian between the qubit and resonator 
can be written as 

Tfq-r = hg fry (a f + a) , (17) 

where g is the coupling strength. 


Supplementary Note 5 Simulation of the Full System Hamiltonian 

Supplementary figure 5a shows a schematic of the configuration for the devices described in this work, where the 
five islands labelled numerically in white text (1-5) follow the electrostatic description of Ref. [1] with island 4 
the center conductor of the coplanar waveguide resonator. Supplementary figure 5b shows an equivalent circuit in 
which islands 1 and 5 have been grounded (equivalent to neglecting the parasitic stripline modes of the resonator). 
Filled black circles in this circuit indicate the five canonical node variables that would be required for its full quan¬ 
tum description. Junction capacitances are labelled Cj a and Cjb for the small and large junctions, respectively, 
and Lg is the geometric inductance of the qubit loop. Supplementary figures 5c and 5d show successive (approxi¬ 
mate) reductions of the circuit to four and three node variables (labelled in blue), respectively. Bold lines in these 
schematics indicate the chosen spanning tree (equivalent to a choice of gauge) [7,8]. 

The quantities in panels (b) and (c) are related to those in (b) according to 


Ct 

C a 


C h 


C sh 


C 13 + 
C 23 + 
C12 + 

C A + 


jr(Ci4 + C15) 


(Cl 4 + C15) 


( 18 ) 


where = C14 + C24 + C34 + C45. The dimensionless capacitive division factor by which the resonator voltage 
is coupled to the qubit is given by: 


Pc — [2C34 (C12 + C25) — 2C24 (C13 + C35)] 

X {2C24C34 + (C24 + ^-X34) [2(C23 + Cja) + Cjb] + (C13 + C I 35)[2(C I 23 + C24 + Cja) + Cjb] 

+ (Cl2 + C25)[2(Ci3 + C23 + C34 + C35 + Cj a ) + Cjb]} • (19) 

Note that the junction capacitances influence this value, in particular when the shunt capacitors are smaller. 

Panel (c) is the minimal circuit that fully captures the experimental qubit devices (excluding the resonator), in 
particular because in nearly all devices the total Csh is determined by the combination of Ct, Cb, and C a . The 
approximation inherent in (d) arises from the fact that the two fundamental Josephson-like modes of the circuit 
(the qubit mode in which nodes 1 and 2 oscillate out of phase, and the “plasma” mode in which they oscillate in 
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phase) couple to the same shunt capacitance C s h- By contrast, in the circuit of panel (c), the plasma mode may 
couple to a larger capacitance Cb + C t > C a than the qubit mode. The result is that, unlike (c) where increasing 
the shunt capacitance leaves the plasma mode frequency far above the region of interest, in the more realistic 
circuit of (d) the plasma mode frequency shifts down as the total shunt capacitance is increased, due to the effect 
of Cb and Ct . For the purpose of the present work, however, it turns out that the circuit of panel (d) gives a good 
approximation (5% or better) to all quantities considered, except inelastic quasiparticle tunneling matrix elements, 
where a correct treatment of tunneling events through the two larger junctions precludes grounding one side of 
those junctions. Note that it is necessary to include the geometric loop inductance (though its effect on the energy 
levels is small) to account for the flux noise matrix elements to this accuracy. 

The Hamiltonians for the two circuits are 


Ti 3 = 4n • E C3 • n + E Ja [l - cos(<p 2 — <Pi + ¥>e)] 

+E Jh [2 - cos((j3 2 - <£ 3 ) - cos(<£i)] (20) 

Ha = 4n • E C 4 • n + E }a [l - cos(<p 2 — <Pi + <Pe)] 

E 

+£jb[2 - COS (if>2 - fi 3 ) - COs(cp 4 - <£i)] + -y (£3 - <^ 4 ) 2 (21) 

where n = {hi} is the vector node charge operator, Ec 3, 4 = e 2 /2 • is the inverse charging energy matrix, 
Ej a and Ejb are the Josephson energies of the small and large junctions, respectively, ip, are the node phase 
operators satifying [tpj. hk] = iSjk , and <p e = 27r < 3> e /$o is the dimensionless external flux through the qubit loop. 
Finally, E L = ($ 0 /2it) 2 / L g is the characteristic inductive energy scale for the geometric loop inductance L g . The 
capacitance matrices for the two circuits are given by: 


C 4 


C 3 


/Cjb + Cja + Csh 

-c Ja - c sh 

V 0 


Cja C s h 

Cjb + Cj a + C s h 
— Cjb 


0 \ 

-Cjb 

Cjb / 


/ Cb + Cjb + Cj a + Ca — Cja — C a 0 

— Cja — Ca Ct + Cjb + Cj a + C a —Cjb 

0 —Cjb Cjb 

V -Cjb 0 0 


-Cjb\ 
0 
0 

Cj b J 


( 22 ) 


(23) 


To diagonalize the full Hamiltonians for the two circuits, we first set the phase across the inductor to zero and 
diagonalize the resulting Hamiltonian (which has either two or three node variables) in a truncated charge basis 
containing the states: —10, —9,..., 1,0,1,..., 9,10 for each island (in units of Cooper pairs), having dimension 

21 2 = 441 or 21 3 = 9261. We use the resulting set of eigenstates and eigenenergies to re-express the 

full Hamiltonian in a product state basis: < 8 > |^) where | v) are linear oscillator states resulting from the loop 

inductance L g and the shunt capacitance across it. The oscillator basis is tmncated at u < 3, and the L g = 0 qubit 
basis is truncated at m < 75 for the three-node circuit, and m < 300 for the four node circuit. The increased basis 
size for the latter case is necessary because the four-node circuit does not have an inductance for its common mode 
(the one in which all islands oscillate together relative to ground), and therefore no potential energy, resulting 
in simple, charge eigenstates. The resulting Hilbert space dimensions for the two circuits are 3 x 75 = 225 
and 3 x 300 = 900, much smaller than that which would have otherwise been required: 21 2 x 3 = 1323 and 

21 3 x 3 = 27783 states. This method is a useful way to efficiently include linear inductances into Josephson 
quantum circuits [9], 

Once the eigenenergies and eigenstates of the qubit circuits are determined, this system is again truncated 
(typically at ~10 qubit energy levels) and then coupled to the resonator using the Hamiltonian 

W q - r = 2e {hi - h 2 ) ■ /3 C Kofa 1 ' + a), (24) 
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where f3 c is obtained from Eq. (19) above, and V r0 = ui r ^'2hZ I is the rms ground state resonator mode voltage for 
a mode impedance Z T . A truncated basis of up to ~ 5 photons in the resonator mode is then used to diagonalize 
the resulting dressed Hamiltonian. 


Supplementary Note 6 High-Power Spectrum and Higher-Level Qubit Tran¬ 
sitions 

To characterize higher levels of the qubit, we perform high-power qubit spectroscopy (Supplementary Figure 6). 

The measurement consists of scanning the drive frequency (cud) as a function of the flux bias with a much higher 
power than that used in the main text in Fig. 2a. Strong driving reveals more qubit transitions in the spectrum, 
including the first four qubit levels, single and multi-photon transitions, and resonator-mediated transitions. The 
same simulation used to develop our noise models is used to add the solid lines in Supplementary Figure 6. To 
match the measured transition frequencies optimally across the entire flux bias range, we generally need to tune 
somewhat the design parameters (Supplementary Table 1), typically by 10—25%. Once adjusted, the fitting is very 
good and reproduces the measured spectroscopy over a wide range of flux and frequency values. 

To characterize the anharmonicity accurately, we use a two-tone low-power pulse technique to measure 0-1 and 
1—2 transitions sequentially (Supplementary Figure 7). A short 7r-pulse at 0-1 transition frequency first prepares the 
qubit at |1) state. Then, a low-power frequency scan resolves the 0-1 and 1—2 transitions. The frequency difference 
gives the anharmonicity which is 490 MHz for Qubit B. This is somewhat less than the 830 MHz predicted solely 
from the design values (Supplementary Table 1). Although the measured value is off by approximately 330 MHz 
from the anharmonicity predicted directly from the design values, it is important to keep in mind that this few- 
hundred MHz difference is on top of qubit frequencies more than 10 times larger (i.e., 4.7-5.2 GHz) and, in this 
sense, the error is rather small. 


Supplementary Note 7 Qubit Parameters 

The parameters for the 22 qubits studied in the work are listed in Supplementary Table 1. 


Supplementary Note 8 Noise Models 

The total T) of the qubit is taken to be: 

1 _ 1 1 1 1 
T\ ~ Tf + jf + + T-f 


(25) 


where the terms come from electric noise, magnetic noise, inelastic quasiparticle fluctuations, and the Purcell 
effect, respectively. 

The coupling of electric and magnetic noise to the qubit is modeled as shown in Supplementary Figure 8, in 
terms of voltage sources SVj weakly coupled to the circuit islands, and a current source SIi weakly coupled to 
the circuit loop. We assume the weak-coupling limit for both of these, where C c i are negligible compared to the 
corresponding node capacitances (diagonal elements of C), and M c is negligible compared to the qubit’s loop 
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inductance L g . In this limit, the electric noise can be expressed as well-defined charge fluctuations SQi for node 
i (with magnitudes independent of the node capacitances), and the magnetic noise as flux fluctuation through the 
loop S&t (with magnitude independent of the loop inductance). The coupling Hamiltonians for these fluctuations 
to the qubit can be written: 


tis® = hS^i (26) 

U SQi = Vi5Qi (27) 


where/; = $o(<^3 — ( P4)/27r/L g is the loop current operator and V; = (Q-C J )j is the voltage operator for node 
i. The resulting contributions to the total decay rate are: 


1 

2 

Tf 

¥ 

1 

2 

t q 

¥ 


|(e|/;|g)| 2 5$(w q ) 

|(e|V|g)| 2 5 Q K) 


(28) 

(29) 


For the Purcell-enhanced decay rate of the qubit excited state, we use the expression: 


= |(ei 7 |(a + a t )|gi 7 )| 2 K 


(30) 


where k is the resonator decay rate and |gi 7 ), |ei 7 ) are the two dressed energy eigenstates of the qubit/resonator 
system in the one-photon subspace. 

As shown in purple in Supplementary Figure 8, the quasiparticle noise contribution to the decay rate is modeled, 
following Ref. 10, as a parallel admittance Y qp and corresponding current fluctuations SI. This approximate 
description is justified since Y qp (u;q) <C Yj(w q ) where Yj is the junction impedance. The resulting decay rate 
is [10,11]: 

i = ^|(e lfc |4e^/ 2 -^t e -^/2|g lfc) |2 5 fc p(Wq) (31) 

1 4 fc= i 

where the sum is over the three junctions, the operator t qp transfers a quasiparticle through junction k, the phase 
offsets are ipk = <Fe for the small junction and zero otherwise, and \gik), \e\k) are qubit energy eigenstates in the 
presence of a single quasiparticle on one electrode of junction k. The effective quasiparticle current noise spectral 
density S qp for junction k is given approximately by [10,11]: 


ck_ 

°qp — x qp ^ 


8Aai 


(32) 


where x qp is the dimensionless quasiparticle density (scaled by the density of superconducting electrons), and 
Aai is the superconducting energy gap of aluminum. The eigenstates of the qubit in the presence of a single 
quasiparticle on each circuit node are obtained using a modified charge representation for that node with 20 basis 
states: —19/2, —17/2,..., —1/2,1/2,..., 17/2,19/2 (in units of Cooper pairs). The Hamiltonian is separately 
diagonalized for each single-quasiparticle configuration using the methods described above. Then the matrix 
elements in Eq. (31) above are evaluated between eigenstates of these different configurations, corresponding to a 
single quasiparticle moving through each of the three Josephson junctions in the loop. The resulting decay rates 
associated with inelastic tunneling through each junction are summed to produce the total rate. Note that we 
neglect processes associated with the presence of two quasiparticles simultaneously, which is justified based on 
our observation of ri qp « 0.26. 
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Supplementary Note 9 Definitions of Power Spectral Densities 

The spectral data shown in Fig. 3b in the main text represent a symmetrized PSD, the Fourier transform of the 
symmetrized autocorrelation function, 

S\(u>) = [ dr exp(-iwr) ^(A(O)A(r) + A(r)A(0)) , (33) 

where A is an operator representing the parameter (e.g., flux, charge, current, voltage, ...) that is fluctuating. 
Clearly, by definition, S\(—uj) = S\(u>), and the effects of noise in the classical (hui q <C k B T) and quantum 
(hoj q k B T) regimes are represented equivalently at both positive and negative frequencies. 

In comparison, the unsymmetrized definition for the PSD, 

/ OO 

dr exp(-iwr) (A(O)A(r)) , (34) 

-OO 

explicitly distinguishes energy absorption and emission by the qubit: noise at negative frequencies corresponds to 
energy absorption by the qubit from its environment, whereas noise at positive frequencies corresponds to energy 
emitted by the qubit to its environment. 

By definition, S\(u>) = \ (.S'/ 1 (w) + S (— w) ), and so either can be used to describe the data presented in this 
work. Here, we elected to use the symmetrized PSD S\(u>), because it has the useful attribute that it can directly 
connect noise data measured in the classical and quantum noise regimes (see Figure 3b in the main text). 

The emission and absorption rates for a system with level-spitting frequency cu q are related to the positive and 
negative part of the spectrum respectively, 


r -=Ej fc A^(cU q ), 

A 

r +=E ( 35 ) 

A 

Here, T_( + ) is the energy emission (absorption) rate to (from) the environment, and k\ is the system noise sen¬ 
sitivity, defined as the derivative of the transversal energy change in units of angular frequency with respect to 
A, i.e., k\ = • It is related to the transition dipole matrix element - as presented in Eq. (2) in the main 

text - by k\ = || (e| D\ \g) . The factor 2 can be understood intuitively as the dipole matrix element address¬ 
ing only one of the two transverse matrix elements, whereas the total transverse energy sensitivity to noise is 
related to both off-diagonal elements. Consider, for example, the traditional flux qubit with the two-level-system 
Hamiltonian H = ftA<x z /2 + hsa x / 2, where he = 2/ p <3>b is the energy bias, I p is the persistent current and 
<f>b is the flux bias [12], While the current operator matrix element | (e| I\ | g) \ « / p , the flux-noise sensitivity 
k$ = | | = 2/p/ft, giving the factor 2. k\ is convenient in practice, as it stands for the change in energy per unit 
change in A. For example, k$ = 2n x 1000 GHz/T>o means e/2n = 1 MHz for a 1 /i<E>o change in flux. 

The decay rate we observe in the inversion-recovery experiment (Fig. 2(c) in the main text) corresponds to the 
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sum of emission and absorption rates, 
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(36) 


This recovers the expression in Eq. (2) in the main text. Note that the relation in Eq. (36) is frequency-independent, 
because it does not differentiate transition direction. To disentangle the up and down rates, one needs the informa¬ 
tion of the equilibrium population or polarization. 

At equilibrium temperature, 


5U (~^) _ cxp f 

S v H) ‘ P V KT) • 


(37) 


This equation indicates that, in the classical, low-frequency limit (hw <C k B T), S u (uj) = S u (—ui) = S(uj). In the 
quantum, high-frequency limit (hco k B T), S v (co) S' u (— co), and so 5 u (w) —> 2 S(co). Therefore, a factor 2 
difference arises at high frequencies between the two PSD definitions. 

Our preference in using symmetrized PSD is due to the fact that Sxfwq) is uniformly related to both exper¬ 
imentally measured dephasing (low-frequency regime) and energy decay T) (high-frequency regime). The 1/f 
noise at low frequencies is classical, arising from an ensemble of fluctuators [13], When extending this 1/f trend 
out to higher frequencies, one does not need to scale by a factor 2 to make the connection. In contrast, when using 
S u (w), a factor 2 is required to make this connection. Intuitively, the origin of the factor 2 arises from the defi¬ 
nition of S u (u>), by which the classical noise power is symmetrically distributed over both positive and negative 
frequencies, whereas the the quantum noise power is entirely captured at positive frequencies. 

Of course, either PSD definition can be made to work, and one must simply be cognizant of the differences 
between the PSD definitions [14] in order to plot the data appropriately. 


Supplementary Note 10 Noise Spectroscopy via Spin Locking 

Supplementary Note 10.1 Spin locking technique 

Spin-locking or T\ p noise spectroscopy is an accurate method developed for resolving noise power spectral den¬ 
sities (PSD) by measuring qubit relaxation rates in the rotating frame during driven evolution. The spectroscopy 
spans the intermediate frequency range, i.e., achievable Rabi frequencies, without substantially undermining the 
locking condition. Details of this method are discussed in Ref. [15], 

When a two-level system (TLS) is driven by a weak (wr <C w q ) and resonant (wr A lu) tone, evolution can 
be conveniently described in the rotating frame, which revolves around the z-axis at the drive frequency (Supple¬ 
mentary Figure 9a). It can be viewed as a fictitious TLS with a quantizing field pointing to X. The level splitting 
is now the (locking) Rabi frequency, wr, rather than w q for the free-evolution case. Note that the corresponding 
longitudinal relaxation time, T\ p , is defined with respect to the new quantization axis. The source of the relaxation 


25 



is noise at the Rabi frequency, transverse to the X-axis. Given the definition in Ref. [15], the rate F lp can be 
expressed as 


Tip = — = -<Sj_x(wr) 

J-lp * 

= y [*Sy(wr) + Sz(wr)] 

= \ [ J-Sx(Wq + Ur) + ^S x (w q “ Ur)] + ^Sz(wr) 

« + |%r) 

= 2^1 + fV ■ (38) 

ffere, Fi = l/T) = |5x(u; q ) is the qubit-frame longitudinal relaxation rate and = IS'z(wr) is the rate 
associated with the Rabi-frequency noise. In Eq. (38), the contributing rotating-frame noise in the second line is 
transformed to qubit-frame noise in the third line. The “l/4”-factor arises from the halved noise amplitude at the 
positive sideband. 

Equation (38) suggests that the noise PSD at the Rabi frequency can be extracted by measuring the qubit-frame 
and rotating-frame longitudinal relaxation rates and making the appropriate subtraction. The spin-locking tech¬ 
nique [16-18] is a straightforward way to measure the rotating-frame longitudinal relaxation. We use the standard 
three-pulse spin-locking sequence (Supplementary Figure 9b), under which the qubit undergoes T\ p relaxation dur¬ 
ing the continuous driving pulse (Supplementary Figure 9c). The recorded decay is fit to an exponential function 
to derive the damping rate Fi p . 

Supplementary Note 10.2 Spectral density for thermal photons in a resonator 

Thermal photons in a resonator with decay rate k have an exponential two-time photon-number autocorrelation 
function C(t) in the small n limit [19,20], 


C(t) = nexp(— kt). (39) 

Correspondingly, by the Weiner-Khinchin theorem, the associated power spectral density (PSD) of the thermal 
photons in the resonator is the fourier transform of Eq. (39), 


SrinM = -j—-2- ( 4 °) 

U) z + K z 

From Bloch-Redfield theory, longitudinal relaxation is connected with the two-time correlation of the transverse 
noise. This also applies to the rotating-frame analogue in the T\ p process. The measured T\ p noise spectrum is 
thus related to the traditional PSD. The effect from the non-Gaussian statistics in the photon noise can be ignored 
in such relaxation process. To derive the effective PSD as seen by the qubit, one must account for the dispersive 
coupling x of the qubit to the resonator and its associated Stark shift (i.e., the frequency shift per photon), and the 
factor rj = ft 2 /(k 2 + 4x 2 ) that scales the effective photon population seen by the qubit. The resulting PSD is 

S, M = (2xf [SnnH}^ m = (2 x) 2 J^ 2 - (41) 
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Supplementary Note 11 Comparing Models for Ohmic Charge Noise and 

Ohmic Flux Noise 

In the main text Fig. 3a (device C), it is clear that our data clearly support the position that 1// flux noise is a T) 
mechanism. However, above 3 GHz in this device, there is an ambiguity between ohmic charge noise (magenta 
dahsed line, Fig. 3a) and ohmic flux noise (grey dashed line, Fig. 3a). We used ohmic charge noise in our models 
in the main text, because there is a known physical basis for its role in relaxation from prior work and it gave a 
slightly better match to experimental results across all 22 devices. In Supplementary Figures 10a - lOf, we apply 
our model to the flux dependence of devices B and C, and to the prediction of T) at the flux insensitive points of 
all 22 qubits under two conditions: 

• using ohmic charge noise as was done in the main text, 

• using ohmic flux noise in place of ohmic charge noise. 

We note that the two models give a similarly reasonable match to the data for all qubits. The one possible exception 
is for device C (Supplementary Figures 10a and 10b) above 6 GHz, where this device is highly sensitive to charge 
noise. In this region, one might make a plausible distinction between the efficacy of the two models. Making a 
stronger distinction between ohmic charge noise and ohmic flux noise will be a topic of future work. 


Supplementary Note 12 Thermal Photon Noise in the Low-Number Regime 

The long-time behavior in Stark shift (Ag 1 J ark ) and dephasing rate (T^ h ) due to thermal photons has a nonlinear 
dependence on the photon number (n) in general. Equations (43)-(44) in Ref. [21] give the dependence for an 
arbitrary ratio between x and «, 


A Stark = | lm[VZ]-^X , 

C = fRe[VZ]-|, (42) 

where Z = (1 + i 2 x/k) 2 + i %x/ K ■ Solving for \/~Z , we have 


Im[ \J~Z] 

R e[VZ] 


— (1 — r 2 ) + \/(l + r 2 ) 2 + 16r 2 n + 16r 2 n 2 


(1 — r 2 ) + ,/(l + r 2 ) 2 + 16r 2 n + 16r 2 n 2 


(43) 


where r = 2 \/k- 

In this work, we are focusing on the situation when n is much smaller than 1. Expanding Eq. (43) to first order 
in n yields: 


Im[ \[Z} = r + —-—^ n , 

1 + 

2r 2 

Re[ \[Z ] = 1 + ——^n . (44) 

1 + r z 
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Substituting these expressions into Eqs. (42) gives Eqs. (4)-(5) in the main text. In fact, the low-number condition 
is r-dependent. The linear approximation is valid when 


n <C 


(1 + r 2 ) 2 
16r 2 


1 

16 


(2 + r 2 



(45) 


Therefore, the condition becomes much looser when the system is in either the strong (r > 1) or weak coupling 
(r < 1) regime, meaning that the linear dependence can extend to higher photon number in both cases. The 
condition is tightest (n<§C 1/4) when r = 1. 


Supplementary Note 13 CPMG pulse sequence and filter functions 

The Carr-Purcell [22]-Meiboom-Gill [23] (CPMG) is a dynamical decoupling pulse sequence that is the multi¬ 
pulse generalization of the Hahn spin-echo [24], The CPMG sequence comprises equally spaced 7r-pulses in 
quadrature (with phases 90°-shifted) with respect to the initial 7r/2-pulse. The technique reduces dephasing due 
to low-frequency noise by a coherent refocusing effect imparted by the tt -pulses. The act of applying 7r-pulses in 
the time-domain can be treated in the frequency domain as a band-pass filter that shapes the noise spectra. The 
passband of this filter is inversely related to the spacing between adjacent 7r-pulses. 

During free evolution, the decay function due to dephasing is written exp[—£(r)], where £(r) is called the 
coherence function and r is the total free-evolution time. Assuming a Gaussian noise environment, the coherence 
function is 

pOO 1 

i{T) = T J 0 (46) 

where S z (lo) is the power spectral density of the longitudinal noise that causes the dephasing, and F(oj,t) is a 
sequence-specific weighting function called the filter function which acts to shape the noise spectrum seen by the 
qubit [25], Assuming infinitely short pulses, the filter function for a CPMG sequence with N (even) 7r-pulses is 
= 4sinc 2 (wr/2) sin 4 (wr/4A^)/cos 2 (u;r/27V). As illustrated in Supplementary Figure 11, this filter 
is essentially a bandbpass filter with a passband that peaks around frequency oj/2tt = N/2r, indicating that more 
7r-pulses will shift the filter to higher frequencies. In addition, the filter bandwidth for a fixed passband becomes 
narrower with larger N. 

Taking into account the effect from finite duration of 7r-pulses, the modified filter function has a general form 
[26,27] 


F^\u>,t) = 


(wr) 


1 + (— l) 1+Ar exp(iwr) + 2 l)"* exp(i uiSjr) cos(a;r 7I -/2) 

i= 1 


(47) 


where Sj e [0,1] is the normalized position of the center of the jth tt -pulse between the two 7r/2-pulses and r n is 
the length of each 7r-pulse, yielding a total sequence length r + Nr n . For 7r-pulses of short duration ( e.g ., t„ = 10 
ns) compared with the total free-evolution time (/.<?., as is typical for our experiment, the bandpass filter 

frequency still peaks near uj/2tt = N/2t. In practice, we use Eq. (47) to find its precise position. 

The photon shot-noise measured in this work has a Lorentzian noise power spectral density centered at zero- 
frequency. The spectral density is essentially frequency independent at low frequencies (the “white-noise” region 
of the Lorenzian), and it decreases at higher frequencies (the “tail” of the Lorenzian). For small N, such that the 
filter passband is in the white-noise region, the dephasing time does not change with N. For large enough N, 
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such that the filter passband reaches the Lorentzian tail region, the noise power contributing to the coherence inte¬ 
gral in Eq. (46) is reduced. For such large N, the dephasing time increases as N increases. Since the transverse 
relaxation time is defined (within a Bloch-Redfield picture) by the rates 1/T 2 = l/2Ti + 1/T^, increasing N in 
this tail region will extend T 2 towards the 2 T\ limit. 

For the qubit described in the main text, T 2j cpmg is approximately 40 fj s for N < 100 and represents a typical 
duration r of the free-evolution. Furthermore, as seen in Fig. 6b in the main text, T^cpmg begins to increase 
around N = 100, reaching T^cpmg = 50 ps at N = 200. Taking r = 50 ps and N = 200, the characteristic 
frequency is 200/(2 x 50 /js) = 2 MHz, consistent with the -3dB point of the Lorentzian spectrum with bandwidth 
n/2n = 1.5 MHz (see Fig. 6a in the main text). For N > 1000, T 2i cpmg saturates at about 85 /us, close to the 
expected value of 2T\. 
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